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Abstract 

We construct energy independent but non-local potentials above inelastic thresholds, in terms 
of Nambu-Bethe-Salpeter wave functions defined in quantum field theories such as QCD. As an 
explicit example, we consider NN — )■ NN + mr scattering processes for n = 0, 1, 2, • • •. We show an 
existence of energy-independent coupled channel potentials with a non-relativistic approximation, 
where momenta of all particles are small compared with their own masses. In the case of two- 
body inelastic scatterings such as AA — )■ AA,A^H,SS, on the other hand, we show that energy- 
independent potentials can be constructed without relying on non-relativistic approximations. We 
also propose a method to extract these potentials using time-dependence of general correlation 
functions. 
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I. INTRODUCTION 



It is important to understand hadronic interactions such as nuclear forces from the point 
of view of their constituents, quarks and gluons, whose dynamics is described by Quantum 
Chromodynamics (QCD). Since the running coupling constant in QCD becomes large at 
hadronic scale, however, non-perturbative methods such as the lattice QCD combined with 
numerical simulations must be employed to investigate this problem systematically. Con- 
ventionally the finite size method [l| has been employed to extract the scattering phase shift, 
but the method is applicable only below the inelastic (particle production) threshold. See 
Refs. 2|, |3| for an extension of this method to multi-channel systems. 

Recently an alternative method was proposed to investigate hadronic interactions and 
it has been successfully employed to extract the potential between nucleons below inelastic 



thresholds! 



ladronic 
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13| and 



sk-Gl. Since then, this method has been applied to other more general 

dn 

interactions such as baryon-baryon interactions ITHI 1 1 . meson-baryon interactions 

n n n □ 

three nucleon forces |14J. 1151 ]. See Refs. [16|, |17[ for reviews of recent activities. 

In the method, called the HALQCD method, a potential between composite particles 
is defined in quantum field theories such as QCD. There are two important properties to 
be proven in quantum field theories, in order to define the potential, which is a quantum 
mechanical object. Let us explain the HALQCD method and these two important properties, 
by considering the A^A^ potential as an example. We first introduce the equal-time Nambu- 
Bethe-Salpeter (NBS) wave function 18|] in the center of mass system defined by 



{0\T {N{r, 0)N{r + x, 0)]} lA^A^, W, Cq); 



(1) 



](0| = in(0| is the QCD vacuum (bra-)state, \NN,W,Co)in is the two-nucleon 



where (0| = 

asymptotic in-state at the total energy W = 2\J + with the nucleon mass niN and 
the relative momentum fc, cq represent quantum numbers other than W such as helicity 
of nucleons and the direction of A;, T represents the time-ordered product, and A^(a;) with 
X = {x,t) is the nucleon operator defined by N{x) = eabc{ua{x)^C'yrjdb{x))qc{x) with the 
charge conjugation matrix C and q{x) = {u{x), d{x))'^. Note that a different choice for N{x) 
is possible as long as N{x) can annihilate the 1-particle nucleon state and the difference 
leads to a difference in NBS wave functions defined from them. Note also that N{x) and 
ip{x) implicitly have spinor and flavor indices. 
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An important property of the NBS wave function for the definition of the potential is 
that, as the distance between two nucleon operator, x = \x\, becomes large, the NBS wave 
function satisfies the free Schrodinger( or equivalently the free Klein-Gordon) equation, 



fc2 _V 

{Ew - Ho) (pw,co{x) ~ 0, Ew = Ho = (2) 

where n = mAr/2 is the reduced mass. In addition, the asymptotic behavior of the NBS 
wave function is described in terms of the phase 6 determined by the unitarity of the S- 
matrix, S = e^*"^, in QCD (or the corresponding quantum field theory). This has been shown 
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originally for the elastic tttt scattering 
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20|, where the partial wave of NBS wave function 



for the orbital angular momentum L becomes 



L . sm{kx- Ltt/2 + 6l{W)) 



- ' " (3) 

as X — > oo at 14^ < W^^ = 4:771^^ (the lowest inelastic threshold). The asymptotic behavior 
of the NBS wave function for the elastic A^A^ scattering has been derived in Ref. 2lJ. The 
asymptotic behavior of the NBS wave function such as eq. (jS]) is the first important property 
that motivates the definition of the potential in QCD. 

The (non-local) potential between two nucleons below the inelastic threshold is defined 
by the equation that 

{Ew - Ho) ipw,co{x) = J d^yU{x,y) (pw,coiy) (4) 

at ly < Wth = 2mfyr + m^. In general, the non-local potential U{x,y) could depend on 
the energy H^jlj]. As we will show, however, an energy-independent potential U{x,y) such 
that eq. (j4]) is satisfied for all W < Wth can be constructed. Therefore, if we solve the 
Schrodinger equation with this potential in the infinite volume, its solutions automatically 
provide correct phase shifts in QCD at all W < Wth by construction. An existence of the 
ly-independent potential U is the second important property to define the potential in the 
HAL-QCD method. 



Using the inner product (/,(?) = / d^x f{x)g{x) where / is the complex conjugate of 
/, we introduce a norm kernel defined by Afwico,W2do = {'^Wi,co,'^W2,do)- Since NBS wave 
functions a.t W < Wth are in general linearly independent^, an inverse A^"^ exists and it 



^ This holds at least for {Wi,co) ^ {W2, do) in the sufficiently large volume. Even if some wave functions 
accidentally become linearly-dependent in small volume, we can remove them, so that our construction of 
the energy-independent potential remains the same. 
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satisfies for Wi,W2 < Wth ^ 

(5) 

W<W,i„co W<Wa„co 

Using tlie inverse norm kernel, we define a ket vector \fw,co) as {x\(pw,co) = fw,co{x) and 
its conjugate bra vector {ipw,co\ as 

SO that the non-local potential can be constructed asjo] 

U= (Ew - Ho)\(pw,co){'^w,co\, (7) 

W<Wti„co 

since it is easy to see that it satisfies the Scrodinger equation (jl]) as 

U\(pw,co) = Y (^Wi ~ ^o) \fWi,do){i^Wudo\fW,co) = {Ew - Ho) \ fw,co) (8) 

as long as W < Wth- It should be noted that the non-local potential which satisfies eq. (jl]) 
at < Wth is not unique. For example, we may add an arbitrary term proportional to 
(1 — P) to the non-local potential U without affecting eq. (jll), where the projection is defined 
P = J2 \^w.co){^w,co\- 

W<Wti„co 

The purpose of this paper is to construct an energy independent (non-local) potential 
which satisfies an appropriate Schrodinger equation at low energy but above inelastic thresh- 
olds in quantum field theories. To make our argument more concrete, we mainly consider 
the A^A^ scattering in this paper. 

In Sec. [ni we demonstrate that energy-independent potentials can be constructed above 
inelastic thresholds if the total energy is small enough such that the non-relativistic approx- 
imation is applicable. In Sec. Ill Al we consider A^A^ — A^A^, NNtt scattering as a simplest 
case, where the total energy W is above 2mjv + ?^7r but below 2mN + 2mT^. In Sec. Ill Bl we 
generalize our construction to larger value of W where the A^A^ — )■ A^A^ + nn scattering for 
higher integer n can occur. In this case, momenta of all particles must be still non-relativistic. 
In Sec. Ill Cl we treat a special case of inelastic scattering such as AB — )■ AB, CD, where non- 
relativistic approximation is not required to construct energy independent coupled channel 

^ We first consider the finite volume, so tliat Wi, W2 take discrete values. We then take the infinite volume 
limit, so that 6wi,W2 and J2w should be replaced by S{Wi — W2} and J dW, respectively. 
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potentials. In Sec. Illlt using results obtained in the previous section, we generalize the time 



dependent method for the extraction of the potential[23| to the case at > Wth, in order to 
treat inelastic processes. Conclusions and discussions are given in Sec. [IVl In Appendix \^ 
we compare the construction of the energy-independent potential above inelastic threshold 
given in the main text with other possible variations. 



II. CONSTRUCTION OF ENERGY-INDEPENDENT POTENTIALS ABOVE IN- 
ELASTIC THRESHOLDS 

We here construct energy-independent (non-local) potentials even above inelastic thresh- 
olds for the A^A^ scattering in the center of mass system. In this report we only consider 
pion productions whose n-th threshold energy is given by = 2mp^ + n x with 
being the pion mass. Extensions to other particle productions such as A^A^ or KK, etc. are 
straightforward. 

We introduce energy intervals defined by A„ = [W^i^, W^"^^) for n = 0, 1, 2, ■ ■ -. Given 
the total energy W , the kinetic energy of the A^A^ + nvr system is denoted by E^^r-, which is 
given by 

where Pi+P2+m=i ki = 0. The corresponding free hamiltonian is denoted by H^. Note that 
E^r cannot be determined form the total energy W alone, except for the elastic scattering at 
n = 0, where E^^ is uniquely determined from a given value of W. Since the determination 
of E^ from W is important to construct potentials from the Schrodingier equation and El^ 
for n > 1 cannot determined from W in general, we restrict our considerations in this paper 
to cases where all momenta Pi,P2, ki, ^2, ■ ■ ■ ,kn are non-relativistic, so that we can write 
W ~ + ioT k = 1,2, ■ ■ ■ ,n at W e A„. (We can exclude k = case since E^r can 
always be determined from W without non-relativistic approximation.) This condition is 
explicitly written as pf < m?^ for i = 1,2 and k1 < for i = 1,2, ■ ■ ■ ,n. Unless otherwise 
stated, we assume this condition in this paper. We roughly estimate how many pions can be 
treated within this approximation. If the total energy of two nucleons with one pion at rest is 
equal to the minimum energy of n-pion production such that Sy^m^ + p'^+uIt^ = 2m]^+nmT^, 
the non-relativistic condition, say p^ ^ 0.9 x m^, leads to n — 1 < ^^(\/7.6 — 2) ^ 5. 
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Therefore we may consider up to A^A^ + Gvr with roughly 5% relativistic corrections. Note 
that some configurations of momenta may become relativistic for a given value of W. We 
exclude such configurations in our consideration of this paper. 



A. Simplest case 

To illustrate our strategy to construct energy-independent potentials, let us consider the 
simplest case at W < W^^ = 2m + 2m^ in this subsection. If G Ai ( 2mAr + <W< 
2m7v + 2m7r ), the inelastic scattering with one pion production (A^A^ — A^A^ + tt) becomes 
possible. We can define in this case a set of 4-independent equal time NBS wave function as 

Zn^w,c,M = (0|r{Ar(a;,0)Ar(a; + a;o,0)}|A^A^,iy,co)i„, (10) 

Z^ZyVJ^^,^(a;o,a?i) = {0\T{N{x,{))N{x + x^,Q)-K{x + xi,{])}\NN,W,co)in, (11) 

Zn^^^^,^{x^) = {Q\T{N{x,Q)N{x + x^,Q)}\NN + Ti,W,ci),,, (12) 

Z^Zy Vi^.ciK, x{) = (0|r{A^(a;, Q)N{x + x^, Q)tx{x + x^, + vr, W, Ci)in, (13) 

where Zjy and Z^^ are renormalization factors for nucleon and pion fields, such that N{x) = 
Zll'^N^{x) and 7r{x) = Zl^'^7r^{x), where N^{x) and 7r'^{x) are renormalized nucleon and 
pion field, respectively. We here consider two asymptotic instates |A^A^, W, co)in and |A^A^ + 
TT, W, Ci)in corresponding to two nucleons and two nucleons plus one pion, where Cq and Ci 
represent quantum numbers other than the total energy W. In the present case, {W, Cq) 
and {W,ci) are equivalent to (si,S2,Pi) and (si, S2, Pi, fci) where Si is the helicity of the 
z-th nucleon and P2 is not independent due to the momentum conservation. As mentioned 
before, W ^ Wq + E^r — Wi + E]y- If distances between all operators become large 



(|cco|, \xi — xq\ — t- oo), we expect (and will indeed show in the separated paper[22]) that 
all NBS wave functions given above satisfy free Schrodinger equations such that 

{e% - H^) ^ 0, (e^ - H^) ^0, z = 0, 1. (14) 

We consider the coupled channel Schrodinger equations for A^A^ and A^A^ + vr, which is 
given by 

«-^o'V^^,c. = J: [ild'ynU'\[x],,[yW^J[y],), e (0,1), (15) 

1=0,1 n=0 
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where [x]o = Xq and [x]i = Xq, Xi. Note that Ely ~ Vl^ — W^^ < ii W e Aq. Our task is 
to show that H^- independent 2x2 potential matrix [/'^^ exists. 

For this purpose, we define vectors from these NBS wave functions at G Ai as 

^w,c.^{^%J[^]o),^l^J[x]i)Y, t = 0,l, (16) 
while at G Aq we take only v^^Vco 

Aco = K,co(Wo),<co(Wl))^ (17) 

where the second component V^vFci([^]i) vanishes as distances between all operators go to 
infinity. (No asymptotic A^A^ + tt state exists at ly < 2m7v + ttIt^.) Note that, instead of 
eq. ( !T7|1 . we may define 

V>°w,co ^ «co(Wo),o)^ (18) 

at G Aq. Since the definition of fw,co at G Aq in eq. (IT7|) will be required in Sec. Illllfor 
the time-dependent method, we use it in the main text of this paper, and the construction 
with eq. f|T8|) and other variations will be discussed in Appendix \M 

As in the elastic case, we introduce the norm kernel in the space spanned by v^^yci 

k 

fc=0,l 1=0 

Here indices i,j run over different ranges depending on values of Wi, W2 such that i G I{Wi) 
and j G I{W2), where I{W) = {0} for W e Aq and I{W) = {0, 1} for W G Ai. Otherwise 
stated, we assume this in this subsection. 

As long linearly independent, the Hermitian operator A/" has an inverse as 

E E i^~')w,c.,We,^^e,,W,,, = ^W.wA.A^ (20) 

Schematically M has a following structure: 

^Aroo(Ao,Ao), Aroo(Ao,Ai), Aroi(Ao,Ai)^ 
M = Aroo(Ai, Ao), Aroo(Ai, Ai), Aroi(Ai, Ai) (21) 
^Ario(Ai,Ao), Ario(Ai,Ai), Arii(Ai,Ai)^ 

where 7V'^^(Aj, Aj) represent a sub-matrix whose components are given by A/^.^^ ^y^,^^ with 
Wi G Aj and PF,- G Aj for i,j,a,b = or 1. The corresponding inverse A/""^ has of course 
the same structure. 
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Using this inverse, we define the ket vector Iv^Vcj) ^^e corresponding bra vector 
{ipwc]y whose k-th components are given by 



{[^WwJ = ^wj[^]k), (22) 



kj 

WiGAoUAi j€l{Wi),dj 



(^kcJWO = E E (A^-')^i^c.^..,¥''v^„.,(W/c) (23) 



for k = 0, 1, where dj runs over states which satisfies non-relativistic condition. It is then 
easy to see that 

k 

k=0,l 1=0 

= ^'^^Wi,W2^ci,dr (24) 
Introducing operators Ew, Hq and U such that 

k 

{[xUiEw - Ho)\[y]i) ^ 6ki{E^ - H^) [] 5^'\^n - 2/J (25) 

n=0 

{[x],\U\[y]i) ^ U'\[x],,\y]il (26) 
the coupled channel Shcrodinger equation (fT5|) can be compactly written as^ 

{Ew-Ho)W\^J = U\^\yJ. (27) 
Now it is easy to construct U which satisfies the above equation as 

u= E E E(^v^-^o)|v'W,c.)(^kcJ, (28) 
vfgAouAi i£i{w) a 



smce 



meAouAi jei(Wi) dj 

An energy-independent potential matrix U indeed exists. Note that U is not unique since, 
for example, one can use eq. (ITS!) instead of eq. (IT7I) for ^y/^di that the resulting potential 
from eq. (!28l) differs from the one with eq. (fT7|) . 

Finally let us consider the Hermiticity of U . A matrix element of U is evaluated as 

UkcM = {^w,jU\^k4) = (Vw.JiEw, - ^o)|^k (30) 



^ Here and hereafter the sum over Ci with i 7^ is alway restricted to non-relativistic states if the number 
of particles is more than 2. 
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while 



{U%,c,,w,<i, = {¥>'w,,,,\{Ew, - Ho)\v\v,J = i^mJiEw, - i^o)|^iv.,d,)- (31) 

Therefore potential U is not Hermite in general. However it is effectively Hermite since in 
practice we solve the Schrodinger equation under the condition that -Ev^i = -£'14^25 which is 
equivalent to Wi = W2 in our non-relativistic approximation. 

B. General cases 

It is not so difficult to extend the argument in the previous subsection to more general 
cases, where the total energy satisfies W < t^^'^m^x+i^ discussed before, the validity of 
the non-relativistic approximation requires rzmax = 5 at most. 

Let us consider ly G Aq U Ai U ■ ■ ■ U A„_^^^^. At G A^ with s < rimax, we define a set 
of the equal time NBS wave functions as 

k 

Z^Z^/'^ip'^Jlx],) = {0\T{N{x,0)N{x + Xo,0)l[7r{x + Xi,0)}\NN + t7r,W,c,),^, i < s, 

1=1 

= 0, i > s, 

(32) 

where indices k,i run from to rimax, but fwcii[^]k) with k > s vanishes, as distances 
among all operators (two nucleons and k pions) becomes large, [x]^ = Xq,Xi, ■ ■ ■ ,Xk and Cj 
represents quantum number other than the total energy W of the instate. In the present 
case, {W, Ci) are equivalent to Si, S2, Pi, ki, ^2, ■ ■ ■ , fcj where si is a helicity of the l-th nucleon. 
The coupled channel Schrodinger equation for this system ed, W E Ag {s < rimax) is given 

by 

(eJ^ - i/o'XJH.) = E / d[y]i u'\[x],, [yWw,M)^ ^ e i{w) (33) 

1=0 

I 

where d[y\i = Y[ d^ym., I{W) = {0,1,- • -,5} for W G A,, and k = 0,1, Note 

m=0 

that c:^ W — W^^ < if A; ^ H^)- is now clear that the non-relativistic condition 
is necessary here to determine E^ from W,Ci if k ^ i. Our task is to show that a W- 
independent (nmax + 1) x (nmax + 1) potential matrix U exists. 



9 



As in the previous subsection, we define vectors of NBS wave functions with (rzmax + 1) 
components at G as 

= Wo), ^]jv,cM^). ■ ■ ■ , y^I^^^r ( Wn..aj)^ , (34) 

where i runs over I{W). 

The norm kernel is defined by 

^w,c,,W2d, = {^m,c^^^W2,dJ = E / c?Wfc<^m,c,(Wfc)<^w2,d,(Wfc), (35) 

k=0 

whose inverse is denoted by J^^'^, where i G I(Wi) and j G I(W2)- The bra and ket vectors, 
defined by 

{[xW^J = <,^(W.), (36) 
(i^kcMk) = E E T.i^-%c.,w.dy^^.4^[^]k), (37) 

Wi j(^I{Wi) dj 

satisfy 

fc=0 

Introducing operators Ev/-, Hq and [/, defined as in eqs. fl25|) and fl26|) . we can construct 

= E E E(^H/-i^o)|^W,c,)(V^kcJ, (39) 
which satisfies the coupled channel equation 

(Ew - Ho)\cp\yJ = U\^\^J. (40) 

It is also easy to see the effective Hermiticity of U: Uy^^^. -^^^. = {U'^ywici,W2dj — ^2 

(with non-relativistic approximation). 

C. Special case without non-relativistic approximation 

In this subsection, we discuss a special case of inelastic scatterings where non-relativistic 
approximation is not required to construct energy independent potentials. Here, coupled 
two-particle scattering channels such as AiBi AjBj with i,j = 1, 2, ■ ■ ■ , rzmax are consid- 
ered. For example, in the baryon scattering in the strangeness S = —2 and isospin 7 = 
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channel, AA, NS and EE appear as asymptotic states if the total energy W in the center 
of mass system is 2mj^ < W < 2mj] + m^r. The method to extract coupled channel poten- 



tials in this kind of situation has already been proposed in Ref. [2J], under an assumption 
that energy independent coupled channel potentials exist. In this subsection we prove this 
assumption. 

Given the total energy W, the relative momentum (squared) and the kinetic energy 
E\^, together with the free Hamiltonian Hq, for AiBi are given by 

,2 V72 



w ^ v^FT^H- v/rfT<, e;,. ^ ^. ff„ ^ ^. K - (41) 

where m^. and m^- are masses of Ai and and m\ is their reduced mass. We here assume 
niAi +mBi < rriAj + for i < j. Note that ifW < IV^i = +^b^, pf and become 
negative. 

We defined NBS wave function for AkBk as 

(Z^,Z5ji/V^,c.(^) = {0\T{Akir,0)Bk{r + x,0)}\AB„W,Ci),^, (42) 

where Z^^,^^^. are renormalization factors defined by Ak{x) = Z^/^AKx) and Bk{x) = 
Z^q'^BI.{x) with bare fields A^, Bk and renormalized fields A^, SJJ, and q represents quantum 
number of the asymptotic instate \AiBi,W,Ci)\ri other than W. The index k always runs 
from 1 to rijnax, while the index i runs over I{W) = 1, 2, ■ ■ ■ , s - 1 if W^^^ <W < H/tl. We 
can show that 

hm (i?^-i7oXc,(^) = 0, (43) 



and Lp^^_[x) carries the information of scattering phase shifts [24 1. 
We define vectors and the corresponding norm kernel as 

{x,k\^\^J = </^^,,^(a;), (44) 

k=i ■' 

where i G I{Wi) and j G I{W2)- Using the inverse M"^ of A/", we construct dual vectors 



which satisfies 



'imax p 



k=i 
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An energy independent nmax x ''^max potential matrix which satisfies the coupled channel 
equation that 

{E!^ - Ho)^'^Jx) = d'yU'\x,y)^^^Jy), (48) 

1=1 •' 

can be constructed as 

U'\x,y) = Yl {E'^-Ho){x,k\cp\yJ{^ljl^JyJ), (49) 
w,i&HW),ci 

which is manifestly energy (W) independent, and is Hermite at fixed W. 



III. TIME DEPENDENT METHOD 



In Ref . 23| , a method to extract hadronic potentials below inelastic thresholds from time 
dependence of correlation functions has been proposed, in order to overcome difficulties in 
the conventional method where NBS wave functions with definite energies are extracted 
from asymptotic behaviors of correlation functions in time. In this section, we extend the 
method so that it can be applicable to the case above inelastic thresholds. 

The normalized correlation function is defined by 

1 ^ 
ZNZ'^/'R\[xUt) = —^{0\T{N{x,t)N{x + Xo,t)Y[ 7l{x + Xi,t)J']yN{0)}\0) (50) 

e '^th* 

for k = 0, 1, 2, ■ ■ ■ ,nmax, where J'nn is some source operator which couples to A^A^ states. 
Inserting the complete set for the A^A^ system that 

1 = E E E 1^^ + Q)i„ i„(A^A^ + 27r, I^, qI + • • • , (51) 

w iei{w) Ci 

where the ellipsis represents states with W > W^'^'' and are neglected hereafter, into the 
above correlation function, we obtain 

R\[x],,t) = E E Ee-^'^^Vw^l^miy^c., (52) 

W i£l{W) Ci 

where 

Aw,c. = in(^^ + c,\JNN{Om, A'M/ = W- Wi ~ E!^. (53) 

Note that R'' automatically contains a sum over W,i & -^(^); which is necessary to define 
the non-local potentials in the previous section but is difficult in practice to perform one by 

12 



one. Note however that states with relativistic momenta may appear in the sum. We here 
assume that contributions from such states can be suppressed by an appropriate choice of 
Jnn- Using the non-relativistic approximation, we can derive 
d 

W,i&I{W),Ci 1=0 

"^max p 

E e-('-^)--*y d[y]iU'\[x],, [y]i) ^ e-^'^VX[?/]0^k 



Ci 

1=0 w,ieHW),Ci 

"^max 



1=0 

We then finally obtain 



"max p 

E e-^'-'^^-' j d[y\iU'^\[x]u. [y\i)R\[y]ut). (54) 



{pi I 'T'max p 

-H^ - -| ■ R\[x],,t) = e'--' E ^-""^7 d[y]iU'\[x\k, [y]{)R\[y]ut), (55) 



23| 



which can be used to obtain combined with the derivative expansion 

We here propose a method to extract U^^ directly. For this purpose, we consider a set of 
more complicated correlation functions defined by 

1 ^ 

R^\[x\uMh-^) = -^777J-m{N{x,t)N{x + Xo,t)X{7l{x + Xrn^) 
e "^'h* m=l 



I 

X /rf3^iV(2/,0)iV(2/ + 2/o,0)n^^(2/ + ys,0)}|0), (56) 

s=l 



which satisfies 

{pi I 'Trnax p 

-i/o' - ^1 ■ R'\Mk. \y]ut) = E e-^^-'^--' j d[z]sU'\[x],, [z]s)W\[z],, [y]utm 

Using real eigenvalues of the Hermitian operator R and their eigenvectors Vm whose /c-th 
component is given by v^{[x\k,t) with m = 0, 1, ■ ■ ■,nmax, we can construct the inverse of 
R as 

^^max 1 

{R~T{MkAy]ut) = E ^€i[^]k,t){iL{[y]i,t)y- (58) 

Note that we remove zero modes with = from R and R~^, so that the dimension of R 
and R~^ are effectively reduced from (nmax + 1) x ('^max + l) to (n^i^x + l—no) x (nmax + l— ''^o) 
where uq is the number of zero modes. 

Using the inverse R~^, we can extract U as 

d 



"■max p I /I I 

u'\[xu [y],) = e-'^-' E I Ms |-^o - Q-A ■ R^'iMk, [z]s,t){R''y\[As, [y]iM 

(59) 
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IV. CONCLUSION AND DISCUSSION 



In this paper, we have shown that energy independent and non-local potentials can be 
constructed from a particular set of NBS wave functions even above inelastic thresholds as 
long as momenta of all particles involved are non-relativistic (Sec. Ill Al and Sec. Ill Bp or 
a number of particles is always two (Sec. Ill Cp . We have also derived a formula to extract 
non-local potentials with non-relativistic approximations using the time dependent method 
proposed in Ref. [23 1. 

By the same method in Sec. Ill Al and Sec. 



II B 



non-local potential for three-nucleon systems 



14 



we can construct an energy independent 



15| and even for systems with more than 



25j . the result in 



3 nucleons. In the case of inelastic scattering such as AA — )■ AA, NE, SS 
Sec. Ill CI has completed the HALQCD method proposed in Ref. [2j], where non-relativistic 
approximation is not required. 

An existence of energy independent potentials, which is one of the important properties 
necessary for the HALQCD method to investigate hadronic interactions, is now established 
in rather general situations. A remaining important property to be proven is an asymptotic 
behavior of NBS wave functions for more than 2 particles and its relation to S'-matrix of the 
corresponding quantum field theory. Results on this issue will be published elsewhere {22!. 
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Appendix A: Comparisons among different constructions 



The energy-independent (rimax + 1) x (rimax + 1) potential matrix in the main text is given 
in the coordinate space by 

t^'' = E E E(^^^-^o)KJ«,J (Ai) 

W i&I{W) Ci 

for < /c,/ < n^ax, where {[xUif'^J = and Efc=o'' W.cjv'^^d,) = 

^^■^Swi,W2^Ci,dy The corresponding coupled channel Schrodingier equation is given by 

{Ew - H,)\cp'^J = Y.U'%wJ (A2) 

1=0 

for < /c < rijnax and i G I(W). 

In this appendix we consider some other constructions of energy-independent potentials 
in terms of NBS wave functions and compare them with eq. (lAip . 

1. Modified wave function vectors 

As already mentioned in the main text, we can define the vectors of NBS wave functions 
using eq. f lTSj) instead of eq. f lT7|) . The corresponding modification to eq. fl32|) becomes 

k 

Z^ZyV^,c,(Wfc) = {0\T{N{x,0)N{x + xo,0)l[7rix + xi,0)}\NN + in,W,Ci)ir., k,t<s, 

1=1 

= 0, otherwise. 

(A3) 

The energy- independent potential U^j (where M represents "modified") is given by the same 
formula in eq. (lAip with modifications by (1A3P to 1^]^^^) ("^vf.cJ; while the corresponding 
Schrodingier equation reads 

iei{W) 



for k, i E liW), where I{W) = 0, 1, 2, ■ ■ ■ , s at G A, = [VT^, 



th 



2. Recursive construction 

We construct another energy-independent potential recursively starting from the potential 
for the elastics A^A^ scattering. 
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We first define the corresponding to A^A^ — )■ A^A^ elastic scattering as 

U'n' = E E(^H^-^o)l¥^°v?,eo)(Ccol (A5) 

where the dual wave function {ipw,co\ \fw,co) satisfies {i'wi,co\¥>W2,do) = ^Wi,W2^co,do at 
Wi, W2 G Aq. This is identical to the elastic potential given in eq. ([7]) and satisfies 

at G Aq. 

We then increase the energy so that W G Ai. A condition that <^wc, < k,i < 1 
satisfy the corresponding Schrdingier equation leads to 

= E E E [(Ew - Ho)\^%J - U'^\<p%j] ii^ll^J, (A7) 

WeAi i=0,l Ci 

where {ipw,ci\ for i = 0,1 satisfy {i^Wi,c,\'^W2,dj) = ^'^^WuW2^c„d, at Wi, W2 G Ai. Note that 
U'^ used here is determined in eq. ( IA6P at the elastic region. We define by imposing 
hermiticity for the potential, i.e. U]^ = {U^^)\ from which we can finally determine 

f^A' = E E E [(Ew - ^o)K,) - U'^Vwj] (i^'wj- (A8) 

weAi 1=0,1 Ci 

We now have at G Aq and 11^ for < 2, j < 1 at G Ai. 

It is not so difficult to extend the above construction to larger W recessively. We assume 
that the s x s potential matrix U^^ is already determined at ly G A^.i for s < n^s.^. At 

G As, f/fcs for /c < s can be obtained by 



u'n = EE 



' W,cJ, (A9) 



/=0 



where (V'vk.cJ for z = 0, 1, ■ ■ ■ , s satisfy {'4^w^,c,Ww24j) = ^'^^Wi,W2Sc„dj at Wi,W2 G A^. 
Using the Hermiticity relation that = {U^^Y for /c = 0, 1, ■ ■ ■ , s — 1, we obtain U'^ as 



Un = EE 

weAs i(^i{w) Ci 



' irwJ- (Aio) 



iEw-HoW^J-Y.U'jl\^wJ 
1=0 

The (s + 1) X (s + 1) potential matrix for < /c, / < s is constructed. We can continue 
this recursive construction until s = n^^x- 

The corresponding Schrodingier equation at G A^ becomes 

{Ew - H,)\^'^J = i:u'rl\v>wj (All) 

for < k,i < s, where off-diagonal elements U^'' ioi k ^ I are Hermite by construction. 
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3. Construction at each energy interval 



We finally give a construction of the potential matrix different at each energy interval. 
At W E As for < s < rimax; the (s + 1) x (s + 1) potential matrix can be constructed as 

f^f = E E E(^^-^o)K,)«j (A12) 

WgAs i€l{W) Ci 

for < A;, / < s, where {i)w,c} for z = 0, 1, ■ ■ ■ , s satisfy Efc=o(^m,cJ<^i^2,rf,) = ^Wi,W2^c,4j 
at VTi, 1^2 e A,. Note that U^^ is identical to f/°° given in eq. (E]). 
The corresponding Schrdingier equation at G becomes 

{Ew - i^o)Ke.) = E UfWw,c.) (A13) 

for < /c, 2 < s. 
4. Comparison 

Properties of the original construction in the main text are as follows. 

1. A size of the potential matrix U'^^ is always (rimax + 1)^ at all W E [W{1^, lyth""""^^)- 

2. A form of the potential matrix given in eq. f lAip is also same at all energy. 

3. We use (rimax + l)-length vectors {Iv^vk.q)' \Vw,c,)^ ' ' ' ^ \¥^]v,cf)}^ which are taken to be 
linearly independent for different values of W e [W{1^, W^'"''^^), i G I(W) and Cj. 

4. The construction can be combined with the time dependent method in Sec. IIIII 
In the case of the modified wave function vectors, we have 

1. A size of the potential matrix t/^j is (s + 1)^ at W E Ag. 

2. The form of ?7|| is same at all energy where ?7|j is defined. 

3. We use vectors {\fwci)y \ fwc,)y ' " " > \^Wc,)y 0, ■ ■ ■ , 0}, which are linearly independent 
for different values of 1^ G [W^^^, PVth""""*"^), i E I{W) and q. The (effective) length of 
these vectors is s + 1 at G A^. 

For the recursive construction, we have 
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1. A size of the potential matrix is (s + 1)^ aXWE A^. 

2. The form of is same at all energy where is defined. 

3. We use l^-Pw^c,)^ which are inearly independent for different values of IV e A^, i G I{W) 
and Cj. 

4. The potential matrix is recursively constructed: At IV e A^, for /c = 0, 1, 2, ■ ■ ■ , s 
are determined from {t/^''^ | k' < k}, while can be obtained from using Her- 
miticity. 

For the construction in Appendix IA3t we have 

1. A size of the potential matrix U^^ is (s + 1)^ at IV G A^. 

2. The form of U^^ is different for each s at IV G A^. 

3. We use the (s + 1) length vectors {\Vwc,)AVwc,)^'''AVwc,)}i which are linearly 
independent for different values of TV G A^, i G /(IV) and q. 

4. Ug^ can be determined at each energy interval A^, without using information of other 
energy intervals. 

We summarize the above properties in table [B 

TABLE I: A comparison of different constructions 



construction 


original 


modified(App. |A l|) 


recursive fApp. lA 2[) 


interval(App. |A 3j) 


size of ?7 at IV G As 


(^max + 1)^ 


(. + 1)2 


(. + 1)2 


{s + lf 


As dependence of U 


no 


no 


no 


yes 


vectors 




{\'Av,c,)}k<s 






feature 


t-dependent method 




recursive 


each As 
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